Azimuthal asymmetry in SIDIS off nuclei as probe for q 
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Abstract 

The energy loss parameter q is one of the fundamental transport parameters of hadronic matter. 
Using the twist-4 collinear approach, we show that the cos (p azimuthal asymmetry in unpolarized 
semi-inclusive deeply inelastic scattering (SIDIS) off a large nucleus at intermediate transverse 
momentum is a sensitive observable for its determination. The effect is due to the suppression of 
the azimuthal asymmetry by final-state multiple scattering. 
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I. INTRODUCTION 



The detailed understanding of the properties of hot and cold nuclear matter, as described 
by its transport parameters is one of the topical problems of QCD, in particular in connection 
with high energy heavy-ion experiments at RHIC and LHC As deeply inelastic scattering 
(DIS) is one of the theoretically cleanest processes in QCD it is tempting to use this probe 
to determine some of them. We argue that this is indeed possible for the transport param- 
eter q if one focusses on transverse momentum dependent observables which are especially 
sensitive to transverse momentum broadening while partons travel through hadronic mat- 
ter. The transverse momentum broadening effect arises from final state multiple parton 
interactions, which are enhanced in nuclear matter. Quite a number of different theoretical 
approaches have been formulated to describe this phenomenon l|-[8| and a variety of precise 
experimental observations will be needed to decide which of these formulations are incorrect 
and which are equivalent. One common parameter appearing in all of these approaches is 
the parton transport parameter q which controls parton energy loss or transverse momen- 
tum broadening squared per unit of propagation length Therefore, the calculation and 
measurement of this transport parameter is an important step toward understanding the 
intrinsic properties of nuclear matter, both cold and hot. 

The multiple parton re-scattering in a large nucleus not only leads to energy loss and 
transverse momentum broadening but also to other nuclear effects, for example, the cos and 
cos 20 azimuthal asymmetries of unpolarized semi-inclusive DIS cross sections. These cos0 
and cos 20 -modulations are sometimes referred to as Cahn effect. In Ref. 9|], it was shown 
that the existence of intrinsic transverse parton momenta in the unpolarized distribution and 
fragmentation functions can generate such modulations at low transverse momentum. Later, 
Cahn effect based descriptions of azimuthal asymmetries were formulated in terms of twist-2 
and twist-3 transverse momentum dependent parton distributions (TMDs) 1CH13|| . At large 
transverse momentum, such asymmetries result primarily from hard gluon radiation 14 1 
which can be calculated in the framework of coUinear factorization. Similar processes are 
also responsible for the azimuthal angle dependence of Drell-Yan dilepton production 15|, 



16|. In the intermediate transverse momentum region, Aqcd P± Q, both, collinear 



factorization and transverse momentum dependent factorization are supposed to apply [Ij 
Indeed, in this special kinematic region, the match between leading power TMD factorization 
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and collinear factorization has been made explicit for the cos 20 asymmetry of the Drell 



Yan lepton pair angular distribution in 



19| . In contrast however, it is known already 



since a while that this equivalence does not extend to sub-leading power TMD factorization 



which suffers from severe problems 
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20| and yields results different from that calculated 



2l|. 



in collinear factorization at intermediate transverse momentum 

Azimuthal asymmetries in SIDIS off nuclei are affected by final state interactions and 
thus provide an alternative way to study properties of the nuclear medium. In fact, the 



nuclear dependence of the angular distribution of Drell- Yan lepton pairs has been ca' 
both in the small x and intermediate x region, using the collinear twist-4 formalism 



culated 



22| and 



the color glass condensate(CGC) model 23|. Nuclear dependent azimutha! 



SIDIS have also been investigated recently based on TMD factorization 
ingredient of the treatment in Ref. 



24 



asymmetries in 



25| . The central 



24l | is the relation between the nucleon twist-3 TMDs and 



the nuclear ones. In this paper, we extend that earlier work to a kinematic region of relatively 
large transverse momenta where the process can be treated within perturbative QCD. To be 
more specific, we use the collinear twist-4 approach to calculate the nuclear dependence of 
the COS0 azimuthal asymmetry in SIDIS at intermediate transverse momentum Aqcd << 
Pj^ << Q , where Pj± is the final state jet transverse momentum, and Q is the virtual 
photon momentum. We restrict ourselves to intermediate transverse momentum and the 
current fragmentation region because here the calculation can be significantly simplified 
with the help of general power counting rules valid in light cone gauge 26|, l27| as we will 
demonstrate in the subsequent section. As result of our explicit calculation, we will show 
that the nuclear dependent part of the asymmetry in this specific kinematic region is directly 
proportional to the parton transverse momentum broadening in a nucleus. 



II. AZIMUTHAL ASYMMETRY IN SIDIS OFF NUCLEUS 



The parton model cross section for the unpolarized semi-inclusive DIS process e(/) + 



p/A{P) — e(/') + J{Pj) + X takes the general form 



12 



29|, 



da 



dxBdzdyd^P^ 



j± 



:i-y + !^)FT + {l-y)FL 



+ (2 - y)^l - ycos0jFcos0j + {l-y) cos(20j)F, 



cos 2(pj 



(1) 
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121 ] ■ We define q = I — I' as the virtual photon 
while Xb = Q^/2P -q , z = P ■ Pj/P ■ q and 



where we use the conventions of Ref. 
momentum and its virtuality as = 
y = P ■ q/ P ■ I are the common DIS variables. The azimuthal angle between the transverse 
momentum of the outgoing parton Pj± and the leptonic plane is denoted by (pj. Four 



structure functions F, depending on Xb, Q'^,z (the fraction of the photon energy carried 
by the jet) and Pj_l, encode the QCD structure of the target and the dynamics of the 
partonic subprocess. It is convenient to use light-cone coordinates for which P^ = P~^p^, 
qi' = -XBP^ + n^Q'^/ilxBP^) withp = (1,0,0,0) and n = (0, 1,0,0). At large Pjx, the four 
functions F can be calculated in collinear twist-2 factorization. As stated above, we restrict 
ourself to the asymmetry at intermediate transverse momentum Aqcd << Pjl << Q in the 
current fragmentation region where power counting rules can be applied. In this kinematic 
region, the power behavior of Ft, Fl, -Fcos^j and Fcos2<pj is l/Pj±, ^/QPjl, 
respectively. The cos0 asymmetry is determined by the ratio between the functions Flj°^'^'' 



and Fuu^Ti which read 
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2l|, 



rptwist—2 



':ptwist—2 
cos (f>j 



1 asCp 



2 In ^5(1 



—1 agCp \ - 2 



p2 

2 In ^5(1 



l + z' 



I 



2z^ 



(2) 



(1-^)+ 

where /i(x) is the normal leading power collinear parton distribution. The index a runs 
over flavors of quarks and antiquarks with fractional charge Cq,. To extract the nuclear effect 
we are interested in, one has to go beyond the leading twist treatment and take into account 
twist-4 contributions. 

The machinery of collinear higher twist factorization was pioneered already in the early 
1980s (SOj, and later frequently applied in hadron spin physics 31| and nuclear physics j3|, 
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32| . The higher twist collinear approach has been well established in both the covariant 



and the light cone gauge. In order to better classify the contributions according to power 
counting rules, we carry our calculation out in the light cone gauge with retarded boundary 
conditions. For retarded boundary conditions, certain collinear twist-4 correlators can be 
directly related to the moment of corresponding TMD distributions. 

Following the standard procedure, the higher twist contributions can be systematically 
isolated by expanding the hard part in the parton intrinsic transverse momentum and in- 
cluding the diagrams with transverse polarized gluon exchange between the struck parton 
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and the target remnant. At twist-4 level, the correlators associated with these two types 
of contributions are of the form {'ipd±d±tp) , {iljd±A±ip) and {tpA±A±ilj) . The general power 



counting rule 26|, |27|] states that the diagrams with one additional transversely polarized 



gluon exchange are suppressed by one additional power of Aqcd/Q in the current frag- 
mentation region where Pj± << Q , as long as final state interactions at = +00 have 



been removed by imposing retarded boundary conditions [28|. Therefore, the possible 
leading powers of the corresponding hard parts convoluted with the correlators {ipd±d±ip) , 
{ipd±A±ip) and {'4'A±A±ip) are suppressed by the powers of ^qcd/ Pj±^ Aqcd/Q^J-L 
^qcd/Q"^ respectively. The explicit calculation shows that the leading power contribution 
from {ipd±d±ip) drops out in the azimuthal angle dependent cross section such that its 
sub-leading part Kq^j-,/ QPj^ generates the non- vanishing cos cf) asymmetry. The hard part 
associated with the correlator {^d±A±il)) contributes to the cos (f) asymmetry with the same 
power Aq^^/QPj^. As a result, we can eventually neglect all diagrams with two trans- 
versely polarized gluon exchanges and are left with the contributions from the correlators of 
the first two types. 

To be more specific, the following three twist-4 coUinear correlators enter the calculation, 

^±ix) = j e^-^^y- {P\mi^{-^d^,)D±Mm\P)d'''' 

^x(x) = j ^ e^-'-^y- {P\mi''l^{-^d±p)DUy)m\P)d''' (3) 

where Dp = —idp + Ap. Since the gauge is completely fixed by choosing retarded bound- 
ary conditions in the light cone gauge, all three correlators can be uniquely brought into a 
gauge invariant form. The nuclear dependence has been encoded in the above twist-4 matrix 
elements. It will become evident when we relate them to the relevant moment of the cor- 
responding nuclear TMDs. The perturbative calculation of the hard coefficients associated 
with these twist-4 correlators is straightforward. The functions Ft and Fcos<j>j can thus be 
expressed as convolutions of the hard coefficients and the twist-4 correlators given above. 
At small transverse momentum Pj^ << Q, the results take a remarkably simple form. 



T^twist-4 _ 1 QsCf ^ 2 



2 In ^5(1 -^) + 



PL ' ' (1-^)- 



:ptwist—A 
cos <f>J 



1 agCp ^ 2 



zPj±Q 27r2 

222(1 + ^2 



21n^5(l-z) 



2^2 



+ 



p2 



Using the QCD equation of motion, one obtains the relation 



- Re(p^{x) + lm^±{x) = xff ^{x) 



where 



dy 

An 



(5) 



(6) 



(7) 



One should note that the calculations just discussed apply to SIDIS of both nuclear and nu- 
cleon targets. All results have the same form and differ only in that the collinear correlators 
are taken inside of a nucleus or a nucleon. The nuclear dependence of the azimuthal angle 
asymmetry is thus generated by the nuclear dependence of the collinear twist-4 correlators, 
which can be best seen from the following relations between the twist-4 correlators and the 
moments of TMDs with retarded boundary conditions. 



= j d-'k^klhixM 
Jf^ix) = I d'k^klhix^h 



(8) 
(9) 

here, fi{x) and f± are the normal leading twist TMD quark distribution and twist-3 TMD 
distribution, respectively. Their matrix element definitions are given by 

dy~d'^y± 



fi{x,k±_) 
f\{.x,ki_) 



2(27r)3 
dy^d'^y^ ■ p+ 



*XBP+y-+*fcxm(p|^(0)7+/:(0,y)V^(?/)|P) 



'''^''^'~+"'^-^^j^{P\mi±P^O,y)^iy)\P) 



(10) 
(11) 



2(27r)3 

For retarded boundary conditions, the transverse gauge links appearing in the above matrix 
elements become unity. The nuclear dependence of the TMD distributions f\{x,kj_) and 



/_l(x, k±) have been worked out and given by 

A 



24 



33| 



ff{x, kj 



71A2F 



dH^e-^^^-'^y/'^^^f^ixJ^) 



A 



kl 



(12) 
(13) 



vrA2F 

where = / d^~q{!^) with the quark energy loss transport coefficient g, which controls 
parton energy loss in a cold nuclear medium and transverse momentum broadening squared 
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per unit of propagation length. The superscripts 'A' and 'N' denote the nuclear and nucleon 
TMDs respectively. Using the relations between nucleon TMDs and nuclear TMDs, we find, 



klffix, k±)d^k± -A J klffix, ki_)Sk^ = 

klf^ix, k^)d'k^ -aJ klf^ix, k^)d'k^ = A/f (a;)A2F 



(14) 
(15) 



where f(^{x) is the ordinary integrated parton distribution function of the nucleon. From the 
above two identities, we can conclude that the difference of the cos (pj azimuthal asymmetries 
is proportional to the amount of transverse momentum broadening. More precisely, when 
z 7^ 1, one has. 



< COS</)j >eA - < COS(j)j >eN-- 



twist— 4 
COS 4>j 



y + yV2 



Tptwist—4 



-^twist—4 
cos lj>J 



eA 



Tptwist—4 



^2F 



(16) 



1-y + 2/2/2 1 + 

This is the main result of this paper, which is valid at intermediate transverse momen- 
tum Aqcd « Pj± « Q- Equation ( fT6l) provides a direct handle to extract the crucial 
parameter q from measurements of the azimuthal asymmetry in nuclei and nucleons. 



III. SUMMARY 



In summary, we calculated the cos (j) azimuthal asymmetry in semi-inclusive DIS off a 
nuclear target within the coUinear twist-4 approach. At intermediate transverse momentum, 
the nuclear dependence of the azimuthal asymmetry is linked to the fc^-moment of the 
two relevant quark TMD distributions. The difference between nucleon TMDs and nuclear 
TMDs is generated by final state interactions, such that the difference in the azimuthal 
asymmetries is sensitive to their strength. To be more specific, the difference between 
the COS0 azimuthal asymmetries in SIDIS off nucleons and nuclei is proportional to the 
transverse momentum broadening in the latter. Therefore, it provides an alternative way to 
pin down the transport parameter q by measuring the nuclear dependence of the asymmetry 
at intermediate transverse momentum. 

The approach we developed in this paper can be extended to study the nuclear dependence 
of azimuthal asymmetries in other processes, such as direct photon production in SlDlS off 
nuclei and Drell-Yan lepton pair production in high energy pA scattering. 
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